Ground-state properties and molecular theory of Curie temperature in the coherent 
potential approximation of diluted magnetic semiconductors 



in 
o 
o 

o 

Q 

o 



-a 
o 

> 
On 
OO 

in 

o 
in 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



Sze-Shiang Fengt, Mogus Mochena 
Physics Department, Florida A & M University, Tallahassee, FL 32307 
(Dated: February 2, 2008) 

Using spin-i description of valence holes and Kondo coupling between local spins and carriers, 
GaAs-based III-V diluted magnetic semiconductors (DMS)are studied in the coherent potential 
approximation(CPA). Our calculated relation of ground-state energy and impurity magnetization 
shows that ferromagnetism is always favorable at low temperatures. For very weak Kondo coupling, 
the density of states (DOS) of the host semiconductor is not modified much. Impurity band can 
be generated at the host band bottom only when Kondo coupling is strong enough. Using Weiss 
molecular theory, we predict a linear relation of Curie temperature with respect to Kondo coupling 
and doping concentration x if the hole density is proportional to x. 
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The ferromagnetism of DMS of III-V-type is not well 
understood. To explain ferromagnetism in DMS, var- 
ious models and approaches have been proposedfl-@- 
Though the models differ from each other in details, they 
all agree that the coupling between the carriers and local 
spins is of fundamental importance. An issue of debate, 
however, is how the exchange between localized spins is 
induced by the carriers. One model for this induced ex- 
change is the Ruderman-Kitttel-Kasuya-Yosida (RKKY) 
interactional 0- Another version which results in con- 
clusion equivalent to RKKY is the Zener modelQ which 
uses the fact that the valence holes are on p-orbitals. A 
third model is the double-exchange (DE) mechanism 0. 
But this model is inconsistent with the charge-transfer 
properties^]. Though RKKY can give a Curie temper- 
ature in agreement with experiment, some argue that 
the RKKY model breaks down herein because the lo- 
cal coupling between the carrier and the impurity spin is 
much larger than the Fermi energy and can not be treated 
perturbatively. In dealing with the effect of the localized 
spins, a key issue is whether or not randomness should be 
taken into consideration. The above models are all mean 
field approximations (MFA) which assume homogeneity 
and neglect randomness. But DMSs are disordered sys- 
tems withpositional disorder of Mn impurities. As con- 
cluded in|2||£|, disorder has a substantial influence upon 
carrier magnetic susceptibility. Hence, any first principle 
consideration should take into account the randomness 
of the impurities. 

A classic method of dealing with randomness is the co- 
herent potential approximation(CPA) 10] which has been 
applied to DMS[n|-[i3. Basing on the formalism of [3 
and assuming very large local spin S while keeping the 
product IS constant (where / is the Kondo -like inter- 
action), obtained the density of states and the re- 
lation between Curie temperature and the doping con- 
centration. Using the averaged carrier Green's function, 
Bouzerar et alll3j arrived at the conclusion that the lo- 



cal coupling between the carrier and the impurity spins 
must be intermediate in order to acquire ferromagnetism. 
In this paper, we use the formalism of CPA in ^5| to 
study the ground-state properties of III-V DMS. In con- 
trast to |12j, we keep S = 5/2 and treat the impurity 
spins fully quantum mechanically. Though it is mostly 
accepted that the effective spin of valence holes is §0- 
|2fj|. we describe here the holes as spin-i fermions. It is 
usually believed that such a description can still catch the 
essential physics. Because of spin-orbit interaction, the p- 
orbitals are spilt (with split-off ~ 0.34eV) into a spin-3/2 
multiplet and a spin- 1/2 multiplet[2l|. Using spherical 
approximation, the kinetic energy of the Luttinger-Kohn 
Hamiltonian 21] for the spin-3/2 multiplet takes the form 
X] Al (S 2 k 2 /2m / _ I )c ] { CA( Ck;^ near the valence top after diago- 
nalization, where m M = m/, ~ 0.5m for /i = ±3/2 and 
— mg ~ 0.07m for fi = ±1/2 (m is the effective 
mass of a free hole ) . The interaction between the spin of 
holes and local 5/2-spins now takes a k-dependent form 
E k ,k' s ' 4 J ( k ' k ')ck' exp(-i(k - k') • R) 18]. Since the 
DOS for parabolic band is g(e) = (l/2vr 2 fi, 3 )(2m) 3 / 2 Ve, 
we have the ratio of DOS gh{s) / gei^) — 19 for heavy 
holes and light holes, i.e., about 95% of valence holes are 
heavy holes. Therefore, it is a valid approximation to 
consider only heavy holes. What is more, since the hole 
density is very small, the Fermi wave vector is supposed 
to be very small and the it is thus a good reasonable ap- 
proximation to consider those k values in the interaction 
term. And this leads to the usual assumption that the 
carriers are shallow holes and the coupling of the shallow 
holes to the Mn 2+ can be described by local Kondo in- 
teraction between 5/2-spins and 1/2-spins. 

We study here one of the most commonly studied 
DMS Gai-jMn^As, where the doping concentration x 
varies from 0.015 to 0.075 in region of interest for 
ferromagnetism 1]. In Ga 1 _ :c Mn a .As, ferromagnetism 
was first realized at Curie temperature of 110K 2J. The 
carriers are holes originating from randomly distributed 
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Mn. The system is highly compensatedj^-ji^l with a 
hole density p only around 10% of the Mn density x. 
There are different kinds of randomness, e.g. substitu- 
tional randomness, interstitial randomness, antisite ran- 
domness, and directional randomness of impurity spin. 
It is commonly agreed now that interstitial Mn atoms 
and antisite As only reduce the hole densities and do 
not affect conduction of holes significantly. Therefore we 
consider only two kinds of randomness, i.e., the random 
substitution of the Mn atoms and the random direction 
of the impurity spins. The model Hamiltonian in our 
description is 

i,j,a i 

where Ui depends on whether i is a Ga or Mn site. For 
Ga-site u\ — uf = Eq J2 a c L c ic > an< ^ f° r Mn-site u; = 
uf 1 = Em J2a c L Cio " + «/icSi • si. Si is the local spin of 
Mn at site i, s = (l/2)cJ.T crcr 'C cr ' is the spin of a hole 
where cJ.(c CT ) is the creation(annihilation) operators for 
holes , spin indices a, a' =|, |, and r = (ti, t%, T3) are the 
three usual Pauli matrices. -Em and Eq are the on-site 
energies for Ga and Mn and are assumed constant. The 
hopping energy iy = t if i,j are nearest neighbors and 
zero otherwise and Jk > is the local Kondo coupling. 
The details of the lattice structure are not crucial in the 
following discussion. According to the general scheme of 
CPA, the virtual unperturbed Hamiltonian is 

^(e) = + (e)) C Lc kCT (2) 

cr,k 

where e is the Fourier frequency variable, is the 

CPA self-energy to be determined self-consistently and 
ik is the Fourier transformation of ty . Then the relative 
perturbation V is given by 

V = H-M'(s) = J2 v i ( 3 ) 

i 

1 



These relations are given in [15J in another context. For 
any f{S z ), the spin average is given by (f{S z )) spin = 

X^ = _s e xs " f(S z )/ Yls*=-S eXS * wriere A is determined 
by the condition (S z ) sp i n = m, m is the given mag- 
netization of the impurity spins. In our single-particle 



where Vi = up = J2a( E G - S CT )c! CT c i(T for Ga and 
«! = v i* = EAEm - S CT )cJ a Ci CT + J K Si • s; for Mn. The 
reference Green's function is (ia\M(e)\ja') = (0\ci a (e — 
^ c ) _1 cj cr , |0) where |0) is the vacuum state of the o 
operators, and the associated i-matrices are ip = i>p/(l— 
@v?),tf = vf/{\ - 3$vf ). So the CPA equation and 
DOS are given by 

(l-z)tp +x(t? I ) spin = (4) 

and 

9*{e) = —ImF^e) (5) 

where (• • -) S pin denotes average over the configura- 
tions of impurity spins and F a (e) = (a\\^\\a) = 
(l/iV)X)k[l/(e - *k - where N is the number of 

lattice sites. As usual, the spin-resolved bare DOS (for 
undoped GaAs) can be approximated by the semicircle 
DOS 

Me) = ^f^ (6) 

where A denotes the half-band width. At zero tem- 
perature, the carrier density for spin a can be ex- 
pressed as n a = j g cr (e)de, where Ef is the Fermi 
energy, and the total carrier density n = + n^. 
The total electronic ground-state energy per site is 
£ 9 = J^ £ [5t( £ ) + 9i(e)]de. Defining V| = E M - 
£ T + {J K /2)S Z ,V L = E M - Ej. - (J K /2)S z ,U a = V a - 
J K /2,W T = (l/4)j2S , -S' + ,VF i = (1/4) J^S+S' , the 
CPA equations can be written as 



(7) 



(8) 

I 

CPA, the Callen-Shtrikman relationj^ that tells there is 
a one-to-one correspondence between m and ((S z ) n ) sp in 
for n > 1 applies . Corresponding to bare DOS, we have 
F<°>(e) = (2/A 2 )(e-Ve 2 - A 2 ) so F a (e) = F(°> (e-£ CT ). 
Solving for £ CT , we have £ CT = e - (A 2 /4)F CT - 1/F a . 



(1-x) 



Eq — £t 
l-F T (£ G -£ 



= 



a-x) 



E G - £4 



1-F 1 (Eq-^ 1 



FIG. 1: DOS and polarization for a number of model param- 
eters 



FIG. 2: Relation of e g versus m. B: J K = 0.1; C:J K = 0.2; 
D: J K = 0.3; E: J K = 0.4; F: J K = 0.5; G: J K = 0.6 



Therefore CPA equations can be turned into equations 
for functions F a for a given e. Once F a (e) are known, 
DOS g&(s) and quantities like Fermi energy can be cal- 
culated. 

To solve CPA equations Q and ©, we chose 
x = 0.05, p = O.lx and set Eq = since we can shift 
the chemical potential without loss of physics. Energies 
are normalized so that A = 1. The solution of F a (e) 
is generated by iteration starting from F 1 -^. First, we 
use the small value of x to generate S CT and then use 
the resulting T, a to generate F a . We calculated the DOS 
and the spin polarization P(e) |2pj | for model parameters 
J K = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, E M = 0,-0.2,-0.4,-0.6 
with 77i = 0.0.5, 1.0, 1.5,2.0.2.5. As a check for our nu- 
merical results, the sum rule J_ g a (z)dz = 1 is pre- 
served and the relation n-p + = p is also preserved 
where n a is calculated from spin-resolved DOS respec- 
tively. In FIG.l we show some of our results. The curves 
presented in the main panels of FIG.l are from the bot- 
tom of the band to Fermi energies at zero temperature 
for clarity since this portion is important for low tem- 
perature physics. The insets outline the full behavior. 
It is seen that for very weak Kondo couplings such as 
Jk = 0.2A,0.3A, there are no impurity bands and the 
DOS is not substantially different from the bare DOS 
in shape. Only when Kondo coupling becomes strong 
enough can there be impurity bands and corresponding 
peaks. In contrast to the conclusion of classical spin 
approximation the spin polarizations are not con- 
stant up to Fermi energies. Fig.l shows that for m > 0, 
there are always more spin-down carriers than spin-up 
carriers, in compliance with the fact that the local p- 
d coupling is antiferromagneticQ. Another substantial 
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FIG. 3: Theoretical estimation of Curie temperature using 
Weiss molecular theory 

difference between our calculated DOS and that in^l is 
that there is no impurity peak above the top of the band 
in our calculation. So the e <-> — e symmetry of the bare 
DOS is broken by the impurities. 

FIG. 2 shows the relation between ground-state energy 
per site and the impurity magnetization m. In actual cal- 
culation of ground-state energy, an important issue is the 
determination of the Fermi energy which is fixed by the 
integration of an interpolation function of the DOS. If 
the interval 8e is not small enough, the integral of the 
interpolated DOS may vary significantly with the choice 
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of 5e. To make it stable enough, we choose Ss = 5 x 1CP 4 
to interpolate DOS linearly. Our results show that for all 
the chosen values of model parameters, the ground-state 
energy per site always decreases, though very slowly, with 
increase of impurity magnetization. Therefore CPA pre- 
dicts that at very low temperatures, ferromagnetism is 
always energetically favorable for all the model parame- 
ters considered. 

In FIG. 3, we show the dependence of Curie tempera- 
ture on the model parameters and the doping concentra- 
tion. Weiss molecular field theory is employed as follows 
to calculate the Curie temperature. Given m, one can 
calculate DOS and then (s z ). So one can establish rela- 
tion (s z ) — (s z )(m). On the other hand , given (s z ), each 
impurity spin feels an effective field Jk(s z ) and thus we 
have m = SB s ((3h) with h = -J K {s z )(m), = l/k B T 
and B$(x) the conventional Brillouin function. For 
very small to, we have (s z ) ~ —Am with A > and 
we have (3h ~ f3J K Am . So B s (0h) ~ (S + l)/3h/3 
and thus the Curie temperature can be estimated by 
k B T c ~ JkS{S + l)A/3 . For small to, letting F T (z) = 
F(z) + i/j(z)m, Fi(z) = F(z) — vb(z)rn where F{z) is the 
paramagnetic solution, we have 

A(J K ,E u ,x,0)~- f delmvb(e) (9) 

where we have ignored the /3-dependence. As our nu- 
merical results (not shown here) indicate, the chemical 
potential is very close to the zero temperature Fermi en- 
ergy in a wide range of temperatures (/3 > 100), showing 
that the Fermi function can be approximated by the zero 
temperature step function. Since the width of full valence 
band of GaAs is about 10— 12 eV (the width of Tg band ~ 
4 eV)H3, here we take A = 6 eV. The left panel in FIG.3 
shows the relation of Tc versus Jk for different values of 
Em- The curves exhibit a linear relation. For Jk > 0.3, 
Tc is almost independent of Em- The right panel shows 
the dependence of Curie temperature on the doping con- 
centration for Em — and various values of Jk- Here we 
still assume that the hole density p = OAx. Again, the 
relation indicated is linear in the range 0.01 < x < 0.07. 
Fig. 3 suggests that to reach the observed Tq — 110K, 
the value of J K needs to be 0.3A - 0.4A. For A = 6 eV, 
this value is much larger than the value 1 eV calculated 
in Il7| . But it is still within the possible range proposed 
in 28] . Like in the dynamical mean field study 29] which 
also used semicircle density of states and bandwidth ~ 
10 eV, the crucial issue in our current CPA study is the 
behaviors of the Curie temperature versus the model pa- 
rameters. The resulting numbers of Curire temperatures 
can be scaled up or down depending on the choice of the 
bandwidth. 

To conclude, we summarize our results here. Using 
CPA and treating the impurity spins fully quantum me- 
chanically, we have calculated the ground-state energies 
of GaAs-based III-V DMS for a wide range of model pa- 



rameters. The results show that ferromagnetism is al- 
ways preferable at low temperatures . Unlike the clas- 
sical treatment of the impurity spins, our approach pre- 
dicts that impurity band can arise only at bottom of 
the band of the host system, showing asymmetry caused 
by doping. With the help of Weiss molecular theory of 
ferromagnetism, we obtained a linear relation of Curie 
temperature with respect to Kondo coupling and doping 
concentration. Our results agrees with that inj^]. As 
is known from experiments [l||22|, Tc increases almost 
linearly with x for x < 0.053 and starts to drop when x 
becomes larger. The contradiction might be reconciled 
by the dependence on x of the exchange integral (usually 
denoted as No/3). It is found in [2^ that the absolute 
value N (3 decreases as x increases, a behavior already 
well known to occur in Coi-^Mn^S. Therefore, if there 
is a relation like Jk = —ax + b with a and b > 0, as sug- 
gested by the data given in [2^ , it can be expected that 
the resulting Curie temperature may start to decrease at 
certain point of doping concentration. 
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